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We establish essentially optimal bounds on the complexity of initial-value problems in the 
randomized and quantum settings. For this purpose we define a sequence of new algorithms 
whose error/cost properties improve from step to step. These algorithms yield new upper com- 
plexity bounds, which differ from known lower bounds by only an arbitrarily small positive 
parameter in the exponent, and a logarithmic factor. In both the randomized and quantum 
settings, initial-value problems turn out to be essentially as difficult as scalar integration. 
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1 Introduction 



The complexity of initial-value problems has been well understood for many years in the de- 
terministic settings, see jl] for a survey of the worst-case and asymptotic results. Much less is 
known about the complexity of these problems if randomized methods or quantum computation 
are allowed. Other problems related to initial- value problems, such as summation or integra- 
tion, have already been well studied in both randomized and quantum settings. A survey of 
complexity results for summation and integration can be found in [3]. 

Randomized and quantum complexity of initial- value problems was first studied in 5 . This 
paper showed that a nontrivial speed-up can be achieved by switching from the worst-case de- 
terministic setting to the randomized or quantum computation. Improved complexity bounds, 
the best known up to now, were next shown in In the Holder class of right-hand side 
functions with r continuous bounded partial derivatives, with the rth derivative being a Holder 
function with exponent p, the e-complexity was shown to be O in the ran- 

domized setting, and O ^(l/e) 1 ^ r+p+1 / 2 ^ on a quantum computer (logarithmic factors are here 
neglected) . 

Unfortunately, these upper bounds do not match known lower bounds on the e-complexity. It 
can be shown J5; that the e-complexity is f2 f ( C r +P+ 1 / 2 ) J in the randomized case, and 
£1 ^(1/ 'g^/C^+p+ln on a quantum computer. Hence, aside from logarithmic factors, there is a 
significant difference between the bounds in the exponent. The gap was essentially filled for a 
restricted class of scalar autonomous problems in 6 . For systems of equations, the discrepancy 
between upper and lower bounds remained an open problem. These results are briefly recalled 
in Section 2. 

In this paper, we show almost optimal upper complexity bounds in the randomized and quantum 
settings for a general initial-value problem. The main results are contained in Theorems 1 and 2. 
We recursively define a sequence A\ , A2 , . . . , Ak of randomized or quantum algorithms (which 
are new for k > 3) for systems of initial-value problems. The performance of Ak improves as 
k grows, which is a crucial component of the complexity analysis. The error and cost bounds 
for Ak are shown in Theorem 1. 

Properties of Ak, with k properly chosen, lead to new upper complexity bounds given in Theo- 
rem 2. These differ from the lower bounds by only an arbitrarily small positive parameter 7 in 
the exponent, and a logarithmic factor. The e-complexity turns out to be O ^(l/e) 1 /( r +p+ 1 / 2 -7)^ 
in the randomized setting, and O ( ( l/e)V( r +P+ 1- 7) J on a quantum computer, up to a logarith- 
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mic factor. Hence, in both the randomized and quantum settings, the complexity of initial-value 
problems is essentially established. Roughly speaking, the problem is as difficult as scalar inte- 
gration P]. 

2 Problem definition and known complexity bounds 

In this section we define the problem to be solved, and briefly introduce basic notions on ran- 
domized and quantum settings used in the rest of the paper. We also recall, for later comparison, 
known results on the complexity of initial- value problems in both settings. 

We consider the solution of a system of ordinary differential equations with initial conditions 

z'(t) = f(z(t)), te[a,b], z(a)= V , (1) 

where / : R d -> R d , z : [a, b] — ► H d and r/ € R rf (f(rj) 7^ 0). Nonautonomous systems 
z'{t) = f(t, z(t)) with / : R^ 1 — > H d can also be written in the form (^Q), by adding one scalar 
equation 



" At) ' 
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z'(t) 




. f«t),z(t)) 



with an additional initial condition u(a) = a. 

We shall study the complexity of problem ((J) in the randomized and quantum settings. We 
assume that the right-hand side function / = [f 1 ,...,f d ] T belongs to the Holder class F r,p 
given as follows. For an integer r > 0, p G (0,1], and positive numbers Dq, D\, . . . , D r , H, 
we define 

F r,P = { / : R d ^ R d | f e C r (H d ), l&fiy)] < A, |<| < r, 

\d r P{y)-d r P{z)\<H\\y-z\\P, y,zGK d , j = 1, 2, . . . ,d}, (2) 

where d 1 f 3 represents any partial derivative of order \i\ of the jth component of /. Here, and 
in the rest of this paper, || • || denotes the maximum norm in R rf . 

Let q = r + p denote the regularity parameter for the class F r ' p . We assume that q > 1, which 
assures that / is always a Lipschitz function. 

We now briefly recall basic notions related to the randomized and quantum settings. For a com- 
plete discussion, in particular of the quantum setting, the reader is referred to 
We wish to compute a bounded function I on [a, b] approximating the solution z. Information 
available about the right-hand side / in the construction of I is provided by a subroutine that 
computes values of a component of / or one of its partial derivatives. In the randomized setting, 
the points at which the values are computed can be selected randomly. On a quantum computer, 
by subroutine calls we mean either applications of a quantum query operator for (a component 
of) /, or evaluations of components of / or its partial derivatives on a classical computer. A de- 
tailed discussion of a quantum query operation can be found in 2 . A mapping A that computes 
I based on the available information is called an algorithm. In both randomized and quantum 
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settings, the approximating function produced by an algorithm is random. 
We now define the error of an algorithm in the two settings. Let (fl, X, P) be a probability 
space. Let {x{\ define the uniform partition of [a, b], so that Xj = a + ih, i = 0,1, ... ,n with 
h = (b — a)/n. The approximation I = l(t) to z = z(t) in [a, b] is constructed in the algorithm A 
based on random approximations = yf to z(xi), i = 0, 1, . . . , n, u) G O,. (We assume that the 
mappings uj i— ► yf are random variables for each / 6 F r ' p .) More specifically, I = l u is given by 

r{t)=^,yf,...,yf){t), (3) 

for t £ [a, b], where ip is a certain mapping defining an approximation over [a, b] from discrete 

approximations yf. 

The error of A at / is defined by 

e»(A,f)= sup \\z(t)-l»(t)\\. (4) 

te[a,b] 

We assume that the mapping u \— * e ul (A, /) is a random variable for each / G i ?r ' p . 

In the randomized setting, the error of ^4 in the class -F r,p is given by the maximal dispersion of 

e UJ (A, /), namely, 

e™ nd {A,F r > p )= sup (Ee"(A,/) 2 ) 1/2 , (5) 

where E is the expectation. (The maximal expected value of e ul {A, f) can be considered as well. 
This only changes constants in the results.) 

In the quantum setting, the error of A in the class F r ' p is given probabilistically by 

e quant (A) F r, P) =e quant (A) F r,p ;(5) = gup M { Q | p{ e " /) > « } < <5 }, (6) 

for a given number (5, < <5 < |. 

The cosi of an algorithm A, cost (A), in the randomized or quantum settings is measured by 
the maximal number (with respect to / £ F r ' p ) of subroutine calls (with a proper meaning of a 
subroutine call in the respective setting) needed to compute an approximation. For a given e > 0, 
the e-complexity of the problem in the respective setting, comp rand ( quant ) (F r ' p , e), is defined to 
be the minimal cost of an algorithm A taken among all A such that e rand (q uant ) (^4 ; i? r .p) < e . 
In the quantum setting, to display the dependence on 5, we shall denote the complexity by 
com^ uant (F r ' p ,e,5). 

We now briefly recall known results on the subject. For scalar autonomous problems (0) (with 
d = 1), matching complexity bounds, up to logarithmic factors, were established in It was 
shown, in a slightly restricted class F r ' p , and for a problem restricted to computing z(b) only, 
that 

. , fn\ 1/(9+1/2) / !\2\ 

comp^F^e) = O U -J flog - J (7) 

and 

, . / n\ 1/(9+1/2) \ 

comp rand (F^,e) =0 1- . (8) 
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In the quantum setting 



/ /i \ V(9+i) / 1 i\ i\ 

comp quant (F r, Pj £> ^ = q M log l og _ + l og log _ (9) 



and, for < <5 < |, 

//1\l/(9+l)\ 

comp quant (F r 'V,5) • (10) 

The proof of these results was based on switching to an equivalent scalar nonlinear equation. 
Such a method cannot be used for systems of equations. The best known complexity bounds 
for a general problem are given in [Hj. It has been shown that 

// n 1/(9+1/3) x \ 

comp rand (F r 'V) =0 - log - (11) 

and 

/ /1\ V(9+l/2) / 1 in \ 

comp^ ant (F r >V,5)=0 (-J (log- + log-J . (12) 

Similarly to the case of scalar equations, lower bounds are given by (JSJ) and ()1U() . Hence, there 

is a significant gap in known estimates of the complexity for systems of equations. 

In the next part of this paper we show that the existing gap can essentially be filled up. 

3 Main results 

In this section we present main results of this paper. The first result, Theorem 1, assures 
the existence of a class of algorithms for solving problem (^Q) in the randomized (RAND) and 
quantum (QUANT) settings, which possess certain error and cost properties. These properties 
will allow us to derive the desired complexity bounds. The algorithms are defined constructively 
in the next section. The next result, Theorem 2, is an immediate consequence of Theorem 1. It 
contains new, and almost optimal, complexity bounds for problem (^Q) in the randomized and 
quantum settings. 

In what follows, L denotes the Lipschitz constant for /. 

Theorem 1 For any fcEN and s = 1,2, ... ,k , there exists an algorithm A s for solving 0) 
with the right-hand side f € F r,p that computes for each n £ N and 5 G (0, ^) an approxima- 
tion I s to the solution z, and has the following properties. 

(i) There exist constants Cf and Cf, depending only on the parameters of the class F r,p , a and b 
(and independent ofn), such that 



sup \\z(t)-l s (t)\\ < 

t£[a,b] 



Cfn as if L(b — a) is arbitrary and n > L(b — a)/ In 2 , 

Cf (6 - a)4 +1 n- Q s if L(b - a) < In 2 andneN, 

(13) 
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where 



f g (2 s -l) + 2 s - 1 -l in RAND 
s \ qa + s-1 in QUANT , 1 ; 



and sup{Cf : b — a < M} < oo for any M > 0. 
TTiis ZioWs f/or n > 5) with probability at least 1 — 5. 

(ii) There exists a constant C|, depending only on the parameters of the class F r ' p , a and b 
(and independent of n), such that for n € N the cost of A s is bounded by 

cost(A s ) < C|rA (p k log n + log ^ , (15) 

where 

a f 2 s -l in RAND 

A = ( s in QUANT . (16) 
and sup{Cf : b — a < M} < oo /or am/ M > 0. 

Proof See Section 5. * 

Let us comment on the conditions in Theorem 1 satisfied by the constants Cf and Cf . The 
algorithm A s will be defined by recursive application of A s -±, . . . , A\ on intervals of length 
going to 0. The properties of C{ and C| assure that constants appearing in the error and cost 
estimates do not grow to infinity as the length of the interval tends to 0. 

The error and cost estimates for algorithm A^ in the randomized and quantum versions, with 
a suitably chosen parameter k, play the main role in establishing the complexity of problem (|T|). 
The crucial point is the convergence from above, as k — > oo, of the sequence {Pk/ctk} to the 
optimal exponent, 

' in RAND 



lim '— = t * T 2 



' in QUANT . 



Note that the results for k = 2 agree with those from [H], since /?2/o2 equals to l/(q + 3) in 
RAND, and l/(q + \) in QUANT. 

The following theorem gives essentially optimal upper complexity bounds for problem lfl|). 

Theorem 2 For any 7 € (0, 1), there exist positive constants C\{"f), 6*2(7) an d £ o(t) (depend- 
ing only on 7, the parameters of the class F r,p , a and b) such that for all e E (0, £0(7)) an d 
5 G (0, i) i/ie e-complexity in the randomized and quantum settings satisfies 

/i\ 1/(9+1/2-7) 

comp rand (F r ^, e) < Ci( 7 ) - (17) 



and 



,'i \ 1/(3+1-7) 1 

comp^ ant (^", e, 5) < C 2 ( 7 ) ( -) log ? . (18) 
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Proof Define 

J [log (l/ 7 + 1)1 in RAND 

\ \2/i\ in QUANT . K > 

Consider the randomized setting. We pass from the probabilistic error to the error given by © 
in a usual way, by selecting a suitable 5. Let K denote a positive (deterministic) upper bound 
on e UJ (Ai c , /), depending only on the parameters of the class F r,p , a and b. (One can see from the 
proof of Theorem 1 that such a bound exists.) Then the randomized error © of A k is bounded 
by 



e rand (A k , F r 'P) 2 < K 2 5 + (Cf 



n 



(for n sufficiently large). Selecting 5 = 3e 2 /(4,K 2 ) and n = f(2Cf /e) 1/afc ] , we get that e™ nd (A k , F r '?) < 
e. The cost of A k is then bounded by 



cost(A fc ) = O U -J log -J 



with the constant in the O-notation depending only on and the parameters of the class F r ' p , 
a and 6. By the definition of A;, we have that 

SOL > +1-2 

Hgiicg 

/1\ 1/(9+1/2-7/2) 1 /n 1/(9+1/2-7) 

cost(^) < d( 7 ) (^-J log - < d(7) f - J 

for suitable constants (^1(7), ^1(7), and e sufficiently small. This proves ifTTI) . 
We handle the quantum setting similarly. We estimate the error e quant (Ak, F r ' p , 5) using l|13|l 
and 1)14(1 . select n to assure that e quant (^4fc, F r ' p , 5) < e, and use the cost bound l|15|) . The choice 
{HJ) of k leads to (HE)). " 



Comparing upper bounds from Theorem 2 with lower bounds given in © and (|1U|) . we see that 
they match up to a small parameter 7 in the exponent. In the quantum setting this holds up to 
a logarithmic factor that depends on 5. 

4 Algorithms in the randomized and quantum settings 

We define algorithms whose properties are the subject of Theorem 1. The definition is recursive. 
In both the randomized and quantum settings, we inductively define a sequence of algorithms 
A±, A2, ■ ■ ■ , A s , . . . , A}, for solving problem ©. The algorithm A^, with a properly chosen index 
k, will be our final algorithm. 

The points {x{} define the uniform partition of [a, b], x% = a + ih, % = 0, 1, . . . , n with h = 
(b — a)/n. We shall call n a basic parameter. To display the dependence on the interval over 
which the algorithm is applied and on the basic parameter used, we shall use the notation 
A s = A s ([a,b],n). An approximation to the solution z computed in A s ([a,b],n) is denoted by 
I s = l s {t), t e [a,b]. 
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We start the induction with s = 1. The algorithm Ai([a, b],n) is defined to be Taylor's algorithm 
in [a, b] with step size h. We set yo = rj. For a given yi (an approximation to z(xi)), we let Zi 
be the solution of the local problem 

= /te(*))> 1 e [^i^i+i] > z i( x i) = Vi > ( 20 ) 

and we define y^+i = ^(xj+i). Here, Z°(i) is the truncated Taylor's expansion of for i G 
given by 

/°(i) = X:(l/j!)^ ) N)(*-^) i - 
i=o 

The approximation to z in j4i([a, 6],n) is defined by 

l 1 (t)=l°(t) for t G (21) 

Suppose that the algorithm A s ([a, b],n) is defined for any [a, b] and n. We now describe how to 
get A s+ i([a, b],n) from A s ([a,b],n). 

We first inductively define a sequence of approximations {yi} in ^4 s+ i([a, 6], n), where yi « -z(xj), 
i = 0, 1, . . . ,n. We set yo = ?7- For a given y^, we apply the algorithm A s ([xi, a^+i], m) to the 
problem (|2U[). with a certain basic parameter m to be chosen later on. The approximating 
function obtained in this way is denoted by If = lf(t), t G [£j,:rj+i]. 

The crucial point, which we now describe, is how we obtain yi+x from and If. We divide the 
interval [xj,Xj+i] into ml subintervals of equal length, with end points 

z l j = Xi+ jh, 

with j = 0,1,... , ml, where h = (xj+i — Xi)/(ml). Here I is another parameter that will be 
chosen in the sequel. (The notation I used for the parameter will not be confused with I s = 
I s (t) denoting the approximating function, since the latter always appears with a superscript. 
The distinction is clear from the context.) 
The solution of ()20|) satisfies the identity 

rnl-1 if 1 

dt. 



Zi{x i+ l) = Zi(Xi) + ]T J f(Zi{t)) 



j=0 

For j = 0, 1, . . . , ml — 1 and i = 0, 1, . . . , n — 1, define the polynomial 



w M = E Tsf^Hit (*}))&/ - n(4)) k ■ (22) 
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Then we can equivalently write 

ml-l ' ij t 1 



Zi(x i+ i) = Zi(xi)+Y^ / Wij(lt(t))dt 

3=0 ti 
-. ml—1 \ 

+ h q+1 ml— Y / ga(u)du (23) 
ml r- i. J 

3=0 o 

ml—1 

+ E / (f(zi(t))-f(im))dt, 



3=0 , 



where function g^ is defined for u € [0, 1] by 

f{lf(z)+uh))- WlJ {ll(z) + uh)) 
9ij(u) = • (24) 

(The value lf(zj +1 ) is meant here as the limit lim lf(t) as t — > Zj +1 from the left.) 

To get the formula for j/j+i, we neglect the last right-hand side term in (|23j) . and we approximate 

the penultimate term. We arrive at the formula 



ml-l 1+ 1 

Vi+l = Vi+ E / w ij( l i(t))dt 
+ h q+l mlKPi(f) . 



(25) 



The vector APj(/) is constructed to approximate the mean of ml integrals 

-i ml — 1 J, 

APi(/)w— E / Wi («)d« (26) 

as follows. 

Approximate the integral of gij by the mid-point rule with N knots Uk, i.e., 

j 1 N-l 

o k=0 

Hence, 

mZ-1 j, mi-lW-l 

- E / E **(«*)• (28) 

The (random) vector APj(/) is computed, in the appropriate setting, by applying the optimal 
randomized or quantum algorithm (with repetitions) to approximate each component of the 
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mean of mlN vectors in the right-hand side of (|28[). For a discussion of algorithms for computing 
the mean, see [2], The approximation error is required to be bounded by £1, 



1 mZ-1 AT-l 
j=0 fe=0 



<ei, (29) 



with probability at least 1 — Si. The parameters iV, si and <5i are to be chosen, and will be 
defined soon. 

An approximation P +1 to z over [a, 6] in the algorithm j4 s+ i([a, b],n) is defined by l s+1 (t) = lf(t) 
for t € Xi + i), and = This completes the definition of A s+ i and, by induction, 

of the sequence of algorithms {A s ([a,b],n)} s >i. 

It remains to define the parameters m, I, N, £\ and Si. They are selected to assure the best 
performance of the algorithms, and result from solving a number of auxiliary optimization 
problems. We do not present here the rather technical analysis leading to the proper selection. 
We give instead the resulting values of the parameters, and prove in the next sections the 
correctness of our choice. To accomplish the definition of algorithm A s+ i we make the following 
choice of parameters: 

, , , n / (n 2 ,^ 1 - 4 ,™ 2 *- 1 ) in RAND 

^ N ) = { („,„-!,„-) in QUANT, (30) 

and we set ei = 1/N. The parameter Si is chosen as a function of 5, the basic parameter n and 
k to be 

| l-d-5)!/^- 1 in RAND 
( 1 - (1 -Sfl n in QUANT. 

It is useful to illustrate with an example the choice of the basic parameter in the recursive 

procedure described above. Consider the randomized case with s = 3. We have that m = n 2 , 

i.e., at each step a new basic parameter is obtained by squaring the current one. The algorithm 

As(\a,b\,ri) is thus defined by n applications of ^([xj, a?i+i], n 2 ) for < i < n — 1. Computing 

A2([xi, Xi+i], n 2 ) requires in turn applications of A\ on subintervals of length (x,i + i — Xi)/n 2 = 

(b — a)/n 3 , with the basic parameter m = (n 2 ) 2 = n 4 . At the final step, the function I s is given 

by Taylor's approximations (polynomials) on subintervals of length (b — a)/(n 3 m) = (b — a)/n 7 . 

(For an arbitrary s, intervals in which I s is a polynomial will be established in the next section.) 

Let us remark that the algorithms Ai and A<i defined above are the well known Taylor algorithm 

and the algorithm analyzed in respectively. For k > 3, algorithms A^ are new. 

5 Performance analysis 

This section is devoted to the proof of Theorem 1. We first show error bounds for the al- 
gorithm A s , neglecting all probabilistic considerations related to the bounds. Such issues are 
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considered in the next subsection, as well as the cost of the algorithm. 



5.1 Error bounds 

We shall prove (|13j) . To derive error bounds, we first discuss the regularity of the approximation 
I s computed in A s ([a,b],n), and the regularity of the function gtj given by (|24|). Note that the 
function I s is in general not continuous. It is a polynomial on each subinterval [c,p) of a uniform 
partition of [a, b] (the last subinterval ending with b is a closed one), with length 

(b-a)/n 2a - 1 in RAND and (b-a)/n s in QUANT. (32) 

To prove this, note that the length of each subinterval in which the approximation I 1 in 
Ai([a,b],n) is a polynomial equals to (b — a)/n. Let the length of each subinterval in which 
I s in A s ([a,b],n) is a polynomial be equal to (b — a)/n^ s . Prom the definition, the function If is 
a polynomial on each subinterval of length 

w , a N f (b-a)/n 1+2 ^ in RAND , OQ , 
(b - a) (nm <ps ) = < ), {' , +(h . ^ TTATl ™ (33) 
v " K ' 1 (6-a)/n i+< P s m QUANT. v ' 

Since functions If define Z s , we get the recurrence relation 

, / 1 + 20, in RAND 

0s+1 -jl + s in QUANT (34) 
with 0i = 1. This gives the desired statement (|32(l. 

Function Z| is thus a polynomial on each interval [zj, see (|23|). (|2*l)l and (|23|). We use this 

property to observe that the following fact concerning gij holds. 

Fact 1 There exist positive constants Mi and M<i, depending only on the parameters of the 
class F r,p , a and b (and independent of i, j, Zj, yj), such that 

\\9ij(u)\\ <Mi, «G[0,1], (35) 

hij(ui) - gij(u 2 )\\ < M 2 \ui - u 2 \, ui,u 2 € [0, 1] , (36) 
and sup{ max{Mi, M 2 } '■ b — a < M } < oo for any M > 0. 

Proof The proof follows from arguments used in the proof of Lemma on p. 828 in 5 . We 
replace h := h, xi := Zj, Xi + ± := Zj +1 , I* := If, y* := lf(zj) and wf := in that Lemma. We 
use the fact that function If is a polynomial in [z l j, Zj +1 ) of the form 

l!(t) = J2(W)zW(z*)(t-zi) k , 
k=0 

where z%j satisfies the equation z[j(t) = f(zij(t)), t € [zj,zj +1 ), with the initial condition 
Zij(z l A = c resulting from the recursive definition of the algorithms. 
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We are ready to prove (|13j) . We proceed by induction. Let s = 1. The desired bound in this 
case is a version of the well known error estimate for Taylor's method. The error of Taylor's 
expansion formula satisfies 

\\zi(t)-i^t)\\ <M(t- Xi y +1 , 

for t 6 [ Xi, Xj+i], where M depends only on the parameters of the class F r,p , a, b and is bounded 
for bounded b — a. First estimating the error — yi\\ in the usual way, and then passing to 

the error over [a, b], we arrive at the bounds 



sup \\z{t) - ^(t)!! < 



e Lh^ e L{b-a) _ i)£-iM/j9 + Mh^ +1 if L(b - a) is arbitrary, 



te[a,b] ~ { 5M(b-a)h q if L(b - a) < In 2 . 

Hence, (j!3j) holds with 

= 2(e L(6 ~ a) - 1)L _1 M(6 - a) ff + M(6 - a)* 1 " 1 , = 5M, and a x = q. 

Suppose by induction that (|T3*|) holds for I s , for any [a, 6] and n. Let e, = — yi. By the 

triangle inequality, 

||e;+i|| < \\z(x i+1 ) - Zi(x i+ i)\\ + \\zi(x i+ x) - yi+\\\. (37) 

From the dependence of the solution on initial conditions, we see that the first term is bounded 
by 

\\z(x i+1 ) - Zi(x i+1 )\\ <e Lh \\ei\\. (38) 
To estimate the second term, we subtract (f25]l from (f23|) to get 



/ ml—l \ 

Zi(x i+l ) - y i+ i = h q+1 ml I — E / gijiu) du - AP»(/) 



J=° o 



mi-1 J+ 

E / (/(*(*))-/(*?(*)))*■ 



mi-1 ^'t 1 

+ E 

i=o 

From the definition of APj(/), and the properties of gij stated in Fact 1, we write the bound 

mi— 1 

r-1 



(39) 



^ E / ^(«)d«-APi(/) 



<CW -1 + ei. (40) 



Here, the term CN^ 1 is the estimate of the mid-point rule error, and e% comes from the ran- 
domized or quantum approximation of the mean of mlN vectors, see ()29j) . (The last bound 
holds with probability at least 1 — 5±.) The constant C depends only on the parameters of the 
class F r,p , a and b, and is bounded for bounded b — a. 

The second right-hand side term in (|39|) will be estimated from the definition of lf(t). 
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Let n > L(b — a)/ In 2. Then L(xi+\ — x\) < In 2, so that the second inequality in (|13f) holds for 
Xj+i], m) (with m = n 2 in RAND and m = n in QUANT). We have 

\\f(zi(t))-f(l!(t))\\ <L\\zi(t)-lt(t)\\ <LC{{x H . 1 -x i )*+ 1 m-«'. (41) 

(The coefficient Cf for AgQxi, aJi+i], m) depends on Xi, Xj+x, but, due to the inductive assumption 
on Cf , it is bounded by a constant depending only on the parameters of the class F r ' p , a and b. 
We denote the bound by the same symbol Cf .) 
Putting together relations (|35|) - (|41|) . we can use (|3T|) to see that 



|e m || < e Lh \\ei\\ + h q+1 ml (CN^ + ei 
+LC{{b - a) q+2 n- {q+2) m- a ° , 



(42) 



for i = 0, 1, . . . , n — 1, where eo = 0. Hence, 

p Lhn _-[._, . . 

IN| < ^ — - ( hq+lml [CN^ 1 + £l J + LCf (6 - a^^rr^W^J , (43) 
for i = 0, 1, . . . , n. 

We now pass to the error in [a, b\. For t € [xj, ajj+i], we have 

\W) - l S+ \t)\\ = \\z(t) - l?(t)\\ < \\z(t) - Zi(t)\\ + \\Zi(t) - If 

< e Lh \\ei\\ + C{h q+l m- a ° . 



(44) 



We now use (|13*|) in ifii)) . Since e Lh — 1 > L/i, we get for arbitrary L(b — a) and n > L(6 — a)/ In 2 
that 

sup \\z(t) - l s+1 (t)\\ < 2 ( e L ( fe - a ) - l) L-\b - a) q {mln)- q (CN- 1 + ex) 
te[a,6] v ' (45) 



^ 2e L(b-a) _ ^ Cf(b-a) q+1 



n -(q+l) m -a s 



If the interval [a, 6] is small enough that L(b — a) < In 2, then exp(L(6 — a)) — 1 < 2L(b — a), 
and we obtain 

sup \\z(t) -l s+1 (t)\\ <4(6-a)"+ 1 (mZn)^'2(CiV- 1 + ei) 

te[a,b] (46) 

+ (4L(6 - a) + 1) Cf (6 - a) q+1 n^ q+ ^m- as . 

Bounds 1)45(1 and 1)46(1 allow us to finish the inductive proof of (|13j) . To close the induction, 
we continue with two cases, randomized and quantum, introducing proper parameters m, I, N 
and e\ given in (|30|) . 

Consider the randomized setting. Introducing the parameters we get 



s+if+Ml J Ci +ln as+1 if L(6 - a) is arbitrary and n > L(b - a)/ In 2 

(47) 



tS5]"*^ l " m ~\Cf +1 (b-a) q+1 n~ a ^ HL{b-a) <ln2 and n E N 



where 

a s+ i = min{ (2 S+1 - l)q + 2 s - 1 , q + 1 + 2a s } . (48) 
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Constants C* +1 and C[ +1 are given by 

C s + l = (b- af (2 [e L{ - b -^ - l) L~\C + 1) + {2e L ^ - l) C{ (b - a)) (49) 

and 

C{ +1 = A{C{ + C + 1). (50) 

Since a s is given by Q14|>. the same holds due to for a s+ i, and the induction is closed in the 
randomized setting. 

In the quantum setting, after introducing the parameters from (|30jl . we arrive at the bound 1)471) . 
where only the exponent a s +i differs from that in the randomized setting. It is now given by 

a s+ i = min{ (s + l)q + s , q + 1 + a s } . (51) 

Constants C{ +1 and C[ +1 remain the same as in the randomized setting, see ()49|) and (|5U|) . 
Since a s+ i agrees with (|14j). the inductive proof of (|TH|) is completed in the quantum setting. 



5.2 Probability and cost 

In the previous subsection we neglected probability issues related to the bounds obtained. The 
cost of computations was also not considered. We now show that (fT3|) holds with probability at 
least 1 — 5, and we prove the cost bound lfT5|) . 
We first show that 

bound (O holds with probability at least (1 - <5 1 )^ n -"\ (52) 
where if>(n, 1) = for all n, and 

s-l 

V/ 4 in RAND 

<=i (53) 

in QUANT 

n — 1 

for s > 2 and n > 2. Indeed, for s = 1 the error bound for Taylor's algorithm holds with 
certainty, that is, ip(n, 1) = . For s = 2, (|13j) holds with probability at least (1 — Si) n , so 
that il>(n, 2) = n, which agrees with ((o^|) . Suppose that (f52*|) and ((53*|l hold for some s > 2. The 
error estimate for l s+1 holds true if the error bounds for If and bounds (|29|) are satisfied for 
i = 0, 1, . . . ,n — 1. This holds with probability at least (1 — J 1 )™ , /'( m > s )+ n . Hence, we have that 

■0(n, s + 1) = nif)(m, s) + n. (54) 

This yields that ip(n,s + 1) satisfies (|53|). which follows after introducing m = n 2 in the ran- 
domized case, and m = n in the quantum case. 
Note that for n > 5, we have 

in RAND 
in QUANT . 




4>(n, s) < 
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This, and the choice of Si given in (|5Tj) . assures that 

(1 _ >i_ 6 

for s = 1, 2, . . . , k, as claimed. 

We now prove the cost estimate (fT5]l by induction. For s = 1, the cost is bounded by ci(r, d)n 
classical evaluations of components of / or its partial derivatives, where c\ (r, d) depends only 
on r and d. Hence, @5)| holds with (3\ = I. 

Suppose that (fT31) is satisfied for A s . The costs of computing successive components that define 
the approximation l s+1 in A s+ i are bounded as follows. The cost of computing If for all i is 
bounded by 

C^nvrfi 3 ^J3k log m + log - 
and the cost of Wij for all i, j by 

C2(r, d)nml, 

where C2(r, d) depends only on r and d. Here, C| can be chosen to depend only on the parameters 
of the class F r ' p , a and b, and is bounded for bounded b — a. 

It remains to take into account the cost of computing APj(/), i = 0, 1, . . . , n — 1. We shall use 
the results on computation of the mean of real numbers, those of Mathe jHj in the randomized 
setting, and of Brassard et al. in the quantum setting. Let re = 2 in RAND, and re = 1 in 
QUANT. To compute APj(/) for i = 0, 1, . . . , n — 1, each satisfying (|29|) . we need at most 

Pn min / mlN, ( — ^ 1 log — 



£iJ J Si 

evaluations of / in RAND, or quantum queries on / in QUANT. The constant P depends only 
on the parameters of the class F r ' p , a and b, and is bounded for bounded b — a. The logarithmic 
factor is related to the number of repetitions needed to increase probability of (|29|) . For the 
basic algorithm for computing the mean, (|29|) holds with probability at least |, and for our 
purposes it must be at least 1 — Si. For a discussion of the role of repetitions and computing 
the median of a number of results, see |2J. 
Putting together these bounds, we get 

cost(Aj + i) < C^nm^ 3 (b^ logm + log + C2(r, d)nml + Pnmin j mlN, ( — \ j log — . (55) 

Note that 

log 1 = log 1 _ (1 _ 1 g)1/w/ifc < c (log \ + k log n 
where c is an absolute constant. 

Finally, we introduce in (|55|> the parameters m, I, N and Ei defined by ()3U[). 
In the randomized case, we arrive at 

cost(^ s+ i) < C 2 s+1 n^+ 1 U k logn + log ^ , (56) 
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where /3 S+1 = max{ 2/? s + 1, 2 S+1 - 1 } (which agrees with (UHJ)), and C 2 S+1 = 2CJ + c 2 (r, d) + Pc. 
In the quantum case, (|ST)|) holds with /3 s+ i = max{ /3 a + 1, s + 1 } (which again agrees with (|TE|0 
and (7 2 s+1 = Cf + c 2 (r, d) + Pc. 

The inductive proof of (|15|) is thus completed. This also ends the proof of Theorem 1. 



6 Final remarks 

We have established, up to arbitrarily small 7 in the exponent and a logarithmic factor, the 
randomized and quantum complexity of initial-value problems. The matching upper bound is 
obtained by defining a sequence of randomized and quantum algorithms for solving the problem. 
The main issues here are the following: 

• the use of integral identity (|2*3*)) satisfied by the solution, 

• proper application of randomized and quantum algorithms for computing the mean of real 
numbers, and 

• the recursive definition of the algorithms. To compute the approximation at the final 
level k, we pass a number of times through all the proceeding levels. 

Let us summarize the results in various settings. Neglecting 7 and the logarithmic factor, we 
have that the complexity is of order (l/e) 1//( - 9+1 ) in the worst-case deterministic setting with 
linear information 0J, and the same holds in the quantum case. In the worst-case deterministic 
setting with standard information, it is of order (1/e) 1 / 9 , see In the randomized setting, we 
stay between these results with complexity of order (l/e) 1 /^ 1 / 2 ) . 
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